High energy gamma rays from old accreting neutron stars 

P. Blasi 1 * 2 

ABSTRACT 

We consider a magnetized neutron star with accretion from a companion 
star or a gas cloud around it, as a possible source of gamma rays with energy 
between 100 MeV and 10 14 — 10 16 eV. The flow of the accreting plasma is 
terminated by a shock at the Alfven surface. Such a shock is the site for the 
acceleration of particles up to energies of ~ 10 15 — 10 17 eV\ gamma photons are 
produced in the inelastic pp collisions between shock-accelerated particles and 
accreting matter. The model is applied to old neutron stars both isolated or in 
binary systems. The gamma ray flux above 100 MeV is not easily detectable, 
but we propose that gamma rays with very high energy could be used by 
Cherenkov experiments as a possible signature of isolated old neutron stars in 
dense clouds in our galaxy. 

Subject headings: acceleration: shock - stars: neutron, old, accreting - radiation: 
gamma 

1. Introduction 

In this paper we study a mechanism for gamma ray production in old magnetized 
neutron stars (NS's) with accretion from a dense gas cloud around it, or from the stellar 
wind of a giant companion. The NS has been supposed to have a surface magnetic field 
Bs ~ 10 9 — 10 10 G, and slow rotation, so that pulsar-like activity can be neglected. 

Old isolated NS's have already been extensively studied in connection with their 
possible X-ray emission by (Zane et al. 1995a, Zane et al. 1995) . In particular (Zane et 
al. 1995) found that inside a distance 16 — 30 pc from the Earth, about 10 old isolated 
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NS's should be present, with a soft X-ray emission in the 0.2 — 2.4 keV band. A statistical 
analysis of their galactic population has been made by (Blaes and Madau 1993) who found 
that in our Galaxy there are ~ 10 9 old isolated NS's, and about 1% of them are in dense 
clouds. 

The accretion of gas onto magnetized compact objects has also been studied by many 
authors (see for example Borner 1980) in order to get a description of the X-ray emission 
by them, and more specifically the interaction of the accreting gas with the magnetic field 
of the NS has been investigated by (Arons and Lea 1976) and (Eisner and Lamb 1977) who 
calculated the shape of the magnetosphere and studied the way for the plasma to reach the 
star surface. 

(Arons and Lea 1976) claimed that a collisionless shock forms at the magnetosphere 
surface. In the present paper we consider such a shock as the site for acceleration of 
particles up to energies ~ 10 6 — 10 8 GeV and calculate the gamma ray luminosity due to pp 
collisions in the accreting gas, in the energy range between 100 MeV and 10 5 — 10 7 GeV. 

Moreover we investigate the possibility that 300 GeV — 1 TeV gamma rays, if 
detectable by Cherenkov experiments, could represent a signature of old isolated NS's in 
our galaxy. This paper is planned as follows: in section 2 accretion is described both on 
isolated NS's and NS's in binaries; in section 3 the diffusive shock acceleration theory is 
applied to calculate the luminosity in the form of high energy particles at the shock around 
the NS magnetosphere; section 4 is devoted to the calculation of the gamma ray flux from 
pp collisions. Discussion and conclusions are presented in section 5. 

2. The accretion 

In this section we consider the accretion of matter onto a magnetized NS in two 
different scenarios: in the first one the NS is isolated and it accretes gas from the interstellar 
medium (ISM); in the second one the NS is in a binary system with a giant star, and the 
accretion works by the stellar wind. In both cases we shall study the role of the magnetic 
field on the accretion flow. 



2.1. Isolated Neutron Stars 

An isolated NS moving with velocity v with respect to the circumstellar medium having 
density p can accrete it with a rate which is given by the Hoyle-Littleton value (Hoyle and 
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Littleton 1939): 



(2GM : 



\2 



M = 71 1 \ I = 3 '° X 10 %) fal (1) 



( V 2 + c |)3/2 

for a flow within the accretion radius: 

2GM S 9 
-R/i = -= y = 1.1 x 10 13 u 5 -o cm 2 

where t>5o = v/(50 km/s) and ni = n/(lcm~ 3 ) and is the sound velocity in the medium. 
A typical value for the velocity of isolated NS's is v ~ 10 — 100 km/s, so that they are 
always highly supersonic with respect to the medium in which they move. Thus in eqs. 
(H) and (fj) we neglected cs compared with v. The radius and the mass of the NS have 
been taken to be Rs = 10 6 cm and Ms = 1 M & respectively. If the sound velocity cs in 
the accreting gas remains appreciably smaller than the free-fall velocity, we can reasonably 
assume, for an ideal gas, that the accretion velocity is 



2GM S 

v ( r )=^—^~ ( 3 ) 

where £ < 1 measures the small possible deviations from the free-fall behaviour, while for 
the density profile it trivially follows that 

holding inside the accretion radius Ra, and in the assumption of quasi spherical accretion, 
with a bow shock at Ra- In this expression for p we took £ = 1 (Bondi and Hoyle 1944) and 
we dropped a factor 4 due to the jump condition at the bow shock usually formed at the 
accretion radius (this shock has been supposed to be a strong one, so that the compression 
ratio is 4). 

Let us now introduce in this accretion scenario the large scale magnetic fields that are: 
i) the dipole-shaped magnetic field produced by the NS 

B NS (r) = B s (^-y s (5) 

where Bs is the surface magnetic field of the star; ii) the magnetic field Bf(r) frozen in the 
accreting plasma. 

During the inflow of the gas towards the NS, the magnetic field lines are bent inside 
the plasma; this bending causes the lines to be compressed according with the conservation 
of the magnetic flux 

B f (r) = B f (R A )(^-r 2 (6) 
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where Ra is given by eq. (|2|). A reasonable assumption is that the magnetic field Bj^Ra) 
at the boundary of the accretion region equals the typical ISM magnetic field ~3x 10~ 6 G. 
The compression of the magnetic field lines predicted by eq. (^) is clearly limited by the 
reconnection rate; more precisely eq. (0) holds up to the point where the hydrodynamical 
time scale becomes comparable with the reconnection time scale: r/v/f ~ t/va, where va 
is the Alfven velocity and f// is the free-fall one. 

The radius at which this happens is given by 

r eq = 7x 10 11 V5o 10/3 nf /3 cm. (7) 

When this radius is reached reconnection compensates the increase of the magnetic field by 
compression and an equipartition value is established for Bf(r). Thus 



B f (R A )(r/R A r 2 r eq <r<R A 



(87T/W 2 ) 1 / 2 r < r e 



The inflow of this magnetized plasma proceeds up to the moment in which the energy 
density (pressure) of the stellar magnetic field equals that of the accreting gas. The radius 
where this condition is fulfilled is usually referred to as the Alfven radius, and in the 
following we shall denote it as Rm, according with the interpretation of this radius as the 
boundary of the NS magnetosphere. 

Two extreme cases are possible: r eq <C R M and r eq 3> Rm- In the first one, at the 
Alfven radius the role of the magnetic field embedded in the accreting plasma can be 
neglected, being much less than the equipartition value. Thus the condition for Rm is 

B% i RM) = p{R M )v\R M ). (9) 

In the second extreme case the plasma reaches the radius Rm already in equipartition with 
Bf(r) and eq. @ becomes 

3dM = 2p{R M )v\R M ). (10) 

The exact equation, taking into account all the contributions is 

B 2 S R M & 2 B 2 f (R M ) 

or("5~) = P\ R m)v {Rm) + — (11) 

57T Ks o7T 

where v and p are given by eqs. (|3]) and (|j) respectively. The physical meaning of the 
Alfven surface and the formation of a shock are widely discussed by (Arons and Lea 1976). 
The problem of describing how the accreting matter is able to reach the surface of the 
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NS is very delicate, and no self-consistent mathematical approach exists for taking into 
account the several processes involved in the interaction of the gas with the magnetosphere. 
Nevertheless some possibilities have been proposed: (Arons and Lea 1976) and (Elsener and 
Lamb 1977) studied the penetration of the gas through the magnetosphere by interchange 
instability. (Borner 1980) and several other authors, mainly in connection with the problem 
of the X-ray emission by NS's, proposed that the plasma accretes down to the surface of 
the NS being channelled along the dipolar magnetic field lines, and forming on the polar 
caps an accretion column. 

In this paper we shall study what happens at the shock, as the site for acceleration of 
protons up to high energies. 

The isolated NS's which we are interested in are those located in dense environments 
(n ~ 10 2 — 10 8 cm -3 ). Following (Blaes and Madau 1993) in our galaxy there should be 
~ 10 9 isolated NS, and ~ 1% of these are inside dense clouds. In the following we shall 
consider the cases n = 100 cm -3 , n = 10 4 cm -3 , and n = 10 8 cm -3 as far as the density of 
the cloud is concerned, and the cases B$ = 10 9 G and B$ = 10 10 G for the surface magnetic 
field of the NS. In table 1 the values of M, Ra, r eq and Rm are reported for Bs = 10 10 G 
and two velocities of the NS: v = 10 km/s and v — 80 km/ s. From these numbers it results 
that the accreting plasma reaches the magnetosphere of the NS already in equipartition 
except in the case of high velocity NS's in very large density clouds (n = 10 8 cm -3 , 
v = 80 km/s). It is easy to see that also in these conditions the equipartition condition 
holds up to v ~ 40 km/s. If equipartition is reached out of the magnetosphere the radius 
Rm comes from eq. fllQl): 

R M = 1.5 x 10 9 v% 7 nf ,7 B% 7 cm (12) 

where B 10 = B S /10 10 G. 

In these calculations we neglected all the effects coming from the rotation of the NS. 
This is equivalent to require that (see Treves et al. 1995): 1) the relativistic wind produced 
beyond the light cylinder is not able to stop the accretion at the radius Ra] 2) the accretion 
flow velocity at Rm is larger than the corotation velocity. 

Condition 1) means, for the period P of the NS, that 

P > 1.1 Bli 2 (M/10 n g/ S y l/4 (R A /W u cm) 1/8 s, 
while condition 2) gives 



P > 4.9(M/10 n g/s)- 3/7 s. 
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It is easy to check that these conditions are fulfilled by old isolated NS's, for which the 
rotation period is usually of a few seconds. 



2.2. Neutron Stars in binaries 

We consider here a specific model for a close binary system with a NS and a giant with 
a stellar wind. The velocity of the wind far from the star is vw = 10 6 — 10 8 cm/s; the rate 
of mass loss will be denoted by Mi oss . The density of wind matter at distance r from the 
giant is 

Pw{r) = (13) 
This radial outflow is appreciably influenced by the presence of the NS at a distance from 
it equal to the accretion radius 

B _^, 7xl0 . (14) 

where Ms, as usual, is the NS mass. We shall use here a simplified model, similar to that 
used in (Berezinsky et al. 1996) (the accreting compact object was a white dwarf there) 
where at ~ Ra a bow shock forms, and inside the shock the accretion onto the NS becomes 
spherically symmetric, with a density given again by eq. (§), where M is now connected to 
Mi oss by geometrical considerations. In particular if d is the interbinary distance, we can 
write 

M ~ M loss A 2 = 7.3 x 10- 8 d^{ y , )~ 4 ( f'T, ) M e /yr, (15) 
V d ' 13K 10 7 cm/s J K 10- 6 M Q /yr J V ' 

where d\3 = d/(10 13 cm). In the following we shall use Mi oss = 10~ 6 M Q /yr. The basic 
features of the accretion and of the interaction between the accreting plasma and the stellar 
magnetic field are the same as those previously explained for isolated NS's. However in 
the case of binaries it is more difficult to fix the boundary conditions on the magnetic field 
at the accretion radius, and some ad hoc assumptions about the stellar magnetic field of 
the companion should be required. For simplicity we shall assume that during accretion 
the equipartition is reached due to the fact that the reconnection time scale becomes 
comparable with the hydrodynamical time scale before the Alfven radius is reached. Thus 
eq. ( [HID holds, and we can write: 

R M = 6.2 x 10 6 B^MZl 17 cm (16) 

where M_ 8 is the accretion rate in units of 10~ 8 M & /yr. 
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In this paper we are interested in NS's that do not show pulsar behaviour, so that, in 
the case of binaries too, we assume that the NS rotates slowly and that its magnetic field is 
not larger than 10 10 G. 

The discussion made in the case of isolated NS's about the formation of a shock as 
a consequence of the interaction of the accreting plasma with the magnetic wall at the 
boundary of the magnetosphere holds in this well; in the case of binaries it is 

possible that the geometry of the accretion is modified with respect to the spherical one, 
with the formation of an accretion disk. This happens when the giant fills its Roche Lobe, 
and the transfer of matter to the NS works by the inner Lagrangian point instead than by 
the stellar wind. In this situation the simplified model used here is no longer valid and the 
effect of the magnetic field should be the disruption of the internal part of the disk itself. 
We shall not consider this case here. 



3. The acceleration 

The shock acceleration mechanism is discussed in several reviews (Jones and Ellison 
1991, Blandford and Eichler 1987, Drury 1983) and we shall not stress here the technical 
details. In this paper we propose that the shock which is formed on the boundary of the 
magnetosphere of a NS, according to the mechanism described in the previous sections, can 
accelerate some fraction of the accreting particles up to very high energies, and that the 
reinteraction of the accelerated particles with the gas produces a gamma ray signal. 



The shock is located at the radius Rm defined by eqs. flT2Tj and (OJ) for the case of 



isolated NS's and NS's in binaries respectively. In the following we shall discuss the two 
cases separately. 

It is easy to estimate the luminosity in the form of accelerated particles (we assume 
they are protons) if we introduce an acceleration efficiency 77: 

gmm 

Lace = V —B • ( 17 ) 

Km 

The value usually used for r\ is ~ 0.1 (Berezinsky et al. 1990) . In the case of isolated NS's, 
by using eqs. ([D and flT2]) we have 

LZ C = 2.6 x 10 26 n v^ 17 B m A/1 nT erg/s. (18) 

For NS's in binaries, by eq. flTE|): 

L Un = L4 x 1Q 37 v ^9/7 ^-4/7 ^ (lg) 
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The maximum energy reachable by this acceleration mechanism can be calculated by the 
condition that the diffusion time of the accelerated particles becomes equal to the typical 
loss time, mainly due, in this case, to inelastic pp collisions, whose typical cross section is 
ct = 3.2 x 1(T 26 cm 2 . 

The diffusion time is given by 

t ml = + Bl) (20) 

Ui - u 2 U\ u 2 

where the subscript 2 refers to the region r < Rm and the subscript 1 refers to r > Rm- 
Di = E (eV) c/ 300 Bi(Gauss) are the diffusion coefficients (we adopt here the Bohm's 
values), and Ui are the flow velocities on the two sides of the shock. In the assumption of a 
strong shock we have U\ju 2 = 4 and eq. (^Cj) becomes 

4 + V2E(eV) R M c 
Uiff- 75 Bs 2GMsR 3 s W 

where we used eq. (0) for B 2 = Bns(Rm) and we used the equipartition condition in the 
accreting plasma. 

The typical loss time is 



, 1 _ 7T^2GM s m p p3/2 

''loss — /n \ — ,V • \^^J 

n(R M )cr c Ma c 
The equation for the maximum energy, obtained equating tdiff and ti oss , is thus: 

E max (eV) = 1.1 x 10 40 ^_^. (23) 

This gives the results shown in table 2 for isolated NS's (except the case shown in the last 
column where the gas reaches the magnetosphere not in equipartition) and 

E b Z x {eV) = 1.8 x 10 15 B w 3/7 Mlf (24) 

for NS's in binaries. 

The differential spectrum of the accelerated particles is taken as N P (E) = K{E + Eq)' 1 , 
where the constant K is easily obtained by the condition 

nax 

dE E N P (E) = L acc , (25) 

where E ~ 1 GeV and E t h is the threshold energy for pion production in pp collisions. The 
exponent 7 in the shock acceleration theory is connected to the compression ratio r = U\/u 2 
by the expression 7 = (r + 2)/(r — 1). For a strong shock r = 4 and 7 = 2. 
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From the previous discussion it results that accreting old NS, both isolated or in 
binary systems, can accelerate particles up to ~ 10 15 — 10 17 eV, by the shock arising on the 
boundary of the NS magnetosphere. The luminosity in the form of accelerated particles 
and reference velocity 50 km/s is shown in table 2. 

4. Gamma ray production 

Gamma rays are produced by pp inelastic collisions between the accelerated particles 
and the accreting gas. The main mechanism for gamma ray production is the decay of the 
neutral pions generated by the reaction p + p — > tt° + X. 

This process has been studied in detail by many authors, and the relevant cross 
sections, taking into account pion multiplicities, can be found in (Dermer 1986) . 

The total number of gamma rays per unit time can be written as follows: 

Q 1 = 2 — / dE N P (E) < fr(E) > (26) 

where x is the grammage out of the magnetosphere and < £o~ > is a multiplicity weighted 
cross section for the reaction p + p — > 7r° + X (Dermer 1986) . 

The real value for x could be appreciably larger than the integral of p onto a straight 
line out to some maximum distance (of the order of the accretion radius), due to the 
presence of the magnetic field, which is responsable for the curvature of the particle 
trajectories and therefore of a larger column density suffered by these particles. 

The integral in eq. (0), if x is simply calculated by integrating over straight radial 
lines, is reported in table 2 for isolated NS's, for several values of the gas density n and 
of the surface magnetic field. For NS's in binary systems (we took M = 10~ 8 M Q /yr, 
d = 10 13 cm) the calculation gives 

Qbin ^ 4 x 1Q 36 £-6/7 g -l ( 27 ) 

It is clear from table 2 that for isolated NS's with v — 50 km/s and for typical distance 
scale of ~ 100 pc, no gamma ray signal can be detected in the region E > 100 MeV for 
n = 100 — 10 4 cm -3 , which is the usual range of average density in nearby dense clouds. 
From this point of view giant molecular clouds (GMC) are of particular interest: the density 
in their cores can reach the value 10 8 cm~ 3 (Turner 1988) so that the gamma ray flux 
above 100 MeV becomes detectable up to a distance ~ 800 pc. In this range of distances 
many giant molecular clouds are present (see Dame et al 1987 ), the closest of which is 
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Taurus-Auriga, at ~ 140 pc. From table 2 we see that an isolated NS in the core of a GMC 
at such a distance should give a signal of 10~ 6 photons cm~ 2 s _1 above 100 MeV. In (Blaes 
and Madau 1993) 19 of such clouds are listed with the correspondent expected number of 
NS's. All of them are inside a distance of 830 pc from the Earth and isolated NS's in their 
possibly very dense cores could produce a detectable gamma ray signal with E > 100 MeV. 

In table 2 it is also reported the flux of gamma rays above 300 GeV that could 
be detected by a Cherenkov imaging experiment such as Whipple, whose sensitivity is 
~ 10~ 13 cmr 2 s~ L . 

While the detectability on scales of 50 — 100 pc from isolated NS's in a cloud with 
n = 10 4 cm -3 is limited to low velocity NS's (v ~ 10 km/s), the signal with E > 300 GeV 
from the very dense cores of the GMC is possible on galactic scales (d ~ 9 — 10 kpc). 



The calculation of the high energy flux has been carried out by introducing the photon 
yield F 7 : 



where x is the grammage and Y 1 ~ 0.1; the photon yield (see Berezinsky et al. 1990) 
represents the number of photons with energy E produced by one proton with the same 
energy which undergoes a pp collision. 



We studied the accretion of matter onto compact objects with surface magnetic field 
B s = 10 10 G and with rotation slow enough to allow accretion and not to show pulsar-like 
activity. 

Such objects could be old NS's, for which it has been foreseen a decrease in the surface 
magnetic field down to 10 9 — 10 10 G (Phinney and Kulkarni 1994) . 

As a result of the interaction between the magnetic field of the star and the accreting 
plasma, a magnetosphere is formed around the NS, bounded by a collisionless shock (Arons 
and Lea 1976), where the pressure of the gas equals that due to the magnetic field of the NS. 
Even if the existence of such a shock has been assumed by several authors, no conclusive 
agreement has been reached about it, due to the fact that a self-consistent mathematical 
description of the accretion in presence of a strong magnetic field does not exist. In this 
paper we assumed the existence of the shock and proposed that it could accelerate nuclei 
(protons) up to high energies (see also Shemi 1995) , according with the usual mechanism 
of the diffusive shock acceleration. 




(28) 



5. 



Discussion and conclusions 
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These accelerated particles can produce a gamma ray signal due to the inelastic 
p + p — > 7T° + X collisions with n° decay in 77, where the target is provided by the accreting 
gas. 

The luminosity in the form of accelerated particles from isolated NS's strongly depends 
on two parameters, the matter density around the NS and the velocity of the NS. The NS's 
we considered here are those contained in dense clouds (there are 10 7 of such objects in 
the Galaxy) where the gas density ranges between 100 cm~ 3 and 10 8 cm~ 3 (in the cores of 
the giant molecular clouds). A statistical study of the velocity distribution of these NS's 
has been made by Blaes and Madau (1993) who found that ~ 6% of these objects have 
v < 20 km/s, ~ 22% have v < 40 km/s and ~ 50% have v < 72 km/s. 

In table 1 we shown the accretion parameters for n = 100, 10 4 , 10 8 cm -3 and for two 
extreme velocities, 10 km/s and 80 km/s, verifying that for old NS's the role of the rotation 
can be neglected. 

In table 2 we shown the acceleration data and the ffuences of gamma rays with 
E > 100 MeV, fixing the velocity at 50 km/s and distinguishing the cases B$ = 10 9 G and 
Bs = 10 10 G, for the usual range of densities. In particular we also estimated the flux of 
very high energy gamma rays (E > 300 GeV) detectable by the Cherenkov experiments 
such as Whipple. From these data it seems hopeless to observe with present detectors the 
100 MeV radiation produced by the pp collisions in gas clouds with n ~ 100 — 10 4 cm -3 , at 
least for distance scale of ~ 100 pc, unless for very low velocity NS's (v < 10 km/s) which 
are however a very small percentage of the total. 

Much more interesting is the situation of the cores of the GMC (n ~ 10 8 cm -3 ): in 
this case the 100 MeV radiation is detectable on distances of 800 pc and the very high 
energy part of the gamma ray spectrum (E > 300 GeV) could be detected by Whipple-like 
experiments on galactic scales (9 — 10 kpc). Many GMC are located within such distances: 
the closest is Taurus-Auriga, at 140 pc, which could contain up to 30 NS's; about 300 NS's 
should be present in Cignus rift, 700 pc far away, and many other GMC can be found within 
800 - 900 pc (Blaes and Madau 1993, Dame et al. 1987) . 

Isolated NS's in the dense cores of GMC would necessarily be bright X ray emitters, 
because the matter stopped at the magnetosphere wall is later channelled to the polar caps, 
liberating a luminosity which is two-four orders of magnitude larger than the gamma ray 
one. However such X rays will hardly leave the dense region, due to the strong Compton 
absorption. Besides accretion of gas down to the star surface could be stopped, for some 
values of the parameters, due to the fact that low velocity NS's in very large density 
environments should accrete at super-Eddington rate: radiation pressure could introduce in 
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this case some non stationary accretion regime. Thus gamma rays would be a unique tool 
for the study of the isolated NS's in overdense regions. On the other hand X rays in the 
0.2 keV region remain the most powerful way to observe isolated NS's in moderately dense 
clouds, where the gamma ray flux should be too small. 

Here we discussed also the case of NS's not showing pulsar-like activity, in binary 
systems, but the conclusions reached about them are strongly parameter and model 
dependent. This is a consequence of two factors: first of all the accretion rate onto the NS 
varies abruptly when the geometry of the binary is changed (e.g. interbinary distance, size 
of the companion), and therefore the position of the magnetic wall is also changed, and 
with it the luminosity in the form of cosmic rays. The second factor is that the geometry of 
the accretion can probably become far from being spherical (disk accretion); in the case of 
disk accretion, the effect of the stellar magnetic field is usually to disrupt the part of the 
disk inside the Alfven radius, and no shock is probably produced in this case, or at least 
it should be not useful for accelerating particles. Nevertheless, in the cases in which our 
calculation holds, the signal from old NS's in binaries can be detected by EGRET if the 
source is located inside ~ 500 pc from the Earth. 

The author thanks V.S. Berezinsky and B.I. Hnatyk for many useful discussions; in 
particular he is grateful to V.S. Berezinsky for having encouraged this calculation. 
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Table 1: The values of the parameters for the accretion onto an isolated neutron star for 
B s = 10 10 G. 





100 


cm 3 


10 4 


cm 3 


10 8 


cm 3 


parameters 


10 km/s 


80 km/s 


10 km/s 


80 km/s 


10 km/s 


80 km/s 


R A (cm) 


2.7 x 10 14 


4.3 x 10 12 


2.7 x 10 14 


4.3 x 10 12 


2.7 x 10 14 


4.3 x 10 12 


M (g/s) 


3.7 x 10 13 


7.3 x 10 10 


3.7 x 10 15 


7.3 x 10 12 


3.7 x 10 19 


7.3 x 10 16 


r eq (cm) 


6.9 x 10 12 


6.7 x 10 9 


3.2 x 10 11 


3.1 x 10 8 


6.9 x 10 8 


6.8 x 10 5 


R M (cm) 


1.0 x 10 8 


6.0 x 10 8 


2.7 x 10 7 


1.6 x 10 8 


2.3 x 10 6 


1.4 x 10 7 



Table 2: The values of the parameters of the acceleration and the gamma ray production by 
pp collisions. 





100 


cm 3 


10 4 


cm 3 


10 8 


cm 3 


parameters 


10 9 G 


10 10 G 


10 9 G 


10 10 G 


10 9 G 


10 10 G 


E m ax (eV) 


3.0 x 10 17 


1.1 x 10 17 


8.1 x 10 16 


3.0 x 10 16 


6.0 x 10 15 


1.6 x 10 15 


L acc (erg/s) 


3.6 x 10 28 


1.0 x 10 28 


1.3 x 10 31 


3.6 x 10 30 


2.0 x 10 36 


5.0 x 10 35 


x (g/cm 2 ) 


1.1 x 10~ 6 


5.6 x 10~ 7 


2.1 x 10^ 4 


1.1 x 10^ 4 


7.8 


4.0 


Q 7 (> 0.1 GeV) (s- 1 ) 


8.0 x 10 21 


1.2 x 10 21 


6.0 x 10 26 


9.0 x 10 25 


7.9 x 10 36 


1.1 x 10 36 


Q 7 (> 0.3 TeV) (s- 1 ) 


8.4 x 10 18 


1.2 x 10 18 


6.2 x 10 23 


9.7 x 10 22 


4.2 x 10 33 


5.8 x 10 32 



